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We review a variety of topics related to hadronic structure functions in exclusive semihadronic tau decays. We 
introduce the concept of structure functions and summarize the most important concepts. We then calculate the 
decay r — > 2>ttv t for very small hadronic invariant mass to one loop in Chiral Perturbation Theory. New interesting 
features emerge with respect to the known results at tree level, in particular for the structure functions Wo and 
we- We discuss the prospects for experimental verification of our predictions. Finally, we discuss various issues 
at higher Q 2 , related to hadronic resonance physics. Here we consider 2n, irK, 3tt, KttK and -kKtt hadronic final 
states. 



1. Overview 

In this note we will present a variety of im- 
portant topics related to differential distributions 
in semihadronic r decays. First, in Sec. 2, we 
will introduce the concept of structure functions 
as a convenient tool to extract the information 
encoded in angular and invariant mass distribu- 
tions. Then in Sec. 3, we will calculate the de- 
cay r — > 3ttis t to one loop in Chiral Perturbation 
Theory. Finally, in Sec. 4, we will discuss various 
issues at higher Q 2 (resonance physics) , consider- 
ing 27r, irK, 37T, KttK and ttKtt final states. 

2. Hadronic Structure Functions 

Consider two and three meson decay modes of 
the t, 



hi{pi)h 2 {p2) 
hi(pi)h 2 (p 2 )h 3 (p 3 ) 



(1) 
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The amplitudes for these decays can be written 

as 

M = Vckm^L^ 
where 



(2) 



and 



(3) 



(4) 



The amplitudes can be parametrized by introduc- 
ing form factors. In the general two meson case, 
there are two such form factors, 

(h 1 h 2 \V 11 - A M |0) = Fv(Q 2 )T^{ Pl ~p 2 ), 



(5) 



In the general three meson case, there are four of 
them. 
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[(pi-p 3 ),F 1 + ( P2 - P3 ) u F 2 }T^ 



j p =i+ 



J p = l- 



(G) 



where the Fi depend on three invariants, e.g. Fi 
Fi(Q 2 ,si,s 2 ). We have used 



LL | LL 

P1+P2 

LL , LL 1 U 

Pi + + 



(P2+P3) 2 and cyclic 



g 2 



(7) 



In specific cases, there are various simplifications. 
If the two mesons are two pions, h\h 2 = 7r~7r°, 
then the vector current is conserved and the scalar 
form factor vanishes, F$ = for m u — m^. 
In the three pion case, h\h 2 h 3 = 7r~7r~7r + or 
7T 7r°7r~, Bosc symmetry relates F 2 to F%, via 
F 2 {Q 2 ,si 7 s 2 ) = Fi(Q 2 1 s 2 , 81). G parity conser- 
vation requires F3 = for m u = ma, and PC AC 
requires F4 = for to„ = n%d = 0. 

So we have seen the the hadronic matrix el- 
ement for two (three) mesons in the final state 
is characterized by two (four) complex functions 
of one (three) hadronic invariants. For experi- 
mental analyses, it turns out to be very useful to 
trade the two (four) complex functions for four 
(sixteen) real valued "structure functions" Wx> 
which are defined from the hadronic tensor 
in the hadronic rest frame jj]], 



:= H^H 1 



(8) 



For the precise definitions, the reader is referred 
to Q, the main points are summarized in Tab. 1. 
The differential decay rate is given by M : 



G 



rfr = iif I^ckm| 2 $>xW0t dPS^ 



(9) 



x 



Note that Wa, Wb and Wsa alone determine 
dT/dQ 2 , 



2\2 



dT{T hihuksVr (M 2 - Q 2 ) 
cx 



dQ 2 



Q 4 



Table 1 

The structure functions 



H v * [ 



J F = 1- 



J = 



J F = 1- 



W A 

WcWdWe 



W F W G 
WhWt 



WsbWsc 
WsdWsg 



J F = r 



W B 



J = 



w 



SA 



hih 2 



hih 2 h 3 



ds\ ds 2 \ (1 



M2 



2 } W^_ 



(10) 



(Almost) all structure functions can be deter- 
mined and disentangled from each other by study- 
ing angular correlations of the hadronic system, 
for details see Q. 

3. Low Energy Expansion of the Decay r — > 

1. Tau decays into three pions involve hadronic 
invariant masses \[Q 2 from threshold 2>M V up 
to M T . Therefore theoretical predictions are dif- 
ficult, and one has to resort to phenomenolog- 
ical models. However, for ^/Q 2 below about 
600 MeV, chiral perturbation theory (CHPT) § 
allows to calculate the hadronic matrix elements 
in a systematic expansion in external momenta 
and quark masses. 

This is interesting for two reasons. First, this 
calculation allows to understand the complicated 
three pion dynamics systematically in a small cor- 
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ner of phase space, which can be used to test and 
improve phenomenological models. Second, the 
small Q 2 region could become accessible experi- 
mentally in future at b and r-charm factories. 

Here we want to emphasize a few new and in- 
teresting features arising from the one loop cal- 
culation in CHPT. We will explain their origin in 
detail, and claim that they should be seen by the 
new high statistics experiments. 

2. Sometimes it is useful to parametrize the 
hadronic matrix element using conventions which 
differ from the notation introduced in Sec. 2. One 
can also use three different form factors F, G and 
H, and introduce isospin indices: 

= (^( Pl )^(p 2 )^(p3)l4(°)l > 

= 6 lj 5 kl 



H 



(pi + p 2 )^G{s 1 ,s 2l s 3 ) 

+ (pi-p2)nH(8l,82,Sa) 

+ p 3 ^F{s\, s 2 , S3) + permutations (11) 

F, G and H are isospin and Lorentz scalars and 
can be written as functions of si, s 2 and S3. 
This general matrix element can be projected on 
a physical matrix element in a trivial manner: 
isospin symmetry allows to relate the matrix el- 
ement of the different decay channels. We have 
calculated them in the framework of CHPT to one 
loop. The complete analytic results and a very 
detailed analysis are given elsewhere ||. Here 
we would like to discuss some interesting features 
of this decay as predicted by CHPT at the level 
of structure functions, showing how they derive 
from basic properties of the form factors. 

We will consider only four structure functions 
Wx, the ones related only to the spin 1 part of 
the matrix element. They can be expressed as: 

( Xl +x 2 ) 2 \K\ 2 



W A = 



2{x\ 



x 2 2 )Re(KH*) 



W c = W A 
W D = Ax 3 (xi 



8x 2 \H\ 2 



x 2 )\H\ 



+ (a?i + x 2 )Re{KH* 
W E = 4x 3 {x 1 +x 2 )lm(KH*) 



(12) 



where simple kinematical functions (see 

Ref. § for details), and K = (G-F)/3. W A gov- 
erns the rate and the distributions in the Dalitz 
plot, while the remaining functions determine the 
angular distribution. 

From CHPT to one loop we find: 



K = 



+ 



i (2 1 
kU + 3FI 

J(s 2 



M 2 , 

s s J(s 3 )-^(J( Sl ) 



1 



4h(a 3 -2M% 



H 



288tt 2 F2 

- 2/2 (si + s 2 - AM I) + \2M% + 2f 6 Q 2 
+ si +s 2 -6(s 3 + M 2 )] +0(/)} 

(5M 2 - 2si)j(si) 



-I — 

5A/ 2 -2s 2 )J(s 2 ) 



f 



96tt 2 F2 



g-2fc)(« 1 



(13) 



where 
J(s) 



16?H 



a In • 



1 



(14) 



a = y/1 - AM 2 /s 

The li are the renormalized coupling constants of 
CHPT at 0(p 4 ) (see Q for numerical values). 

As required by Bose symmetry K (H) is sym- 
metric (antisymmetric) under the exchange of si 
and s 2 . Isospin symmetry requires that the form 
factors in the 27r°7r~ mode are equal to the isospin 
symmetric form factors ( p"3| ) multiplied by a fac- 
tor v2) whereas the form factors in the 2ir~ir + 
mode are given by: 

= K(s 1 ,s 3 ,s 2 ) +K(s 2 ,s 3 ,si) 

- H(si,s 3 ,s 2 ) - H(s 2 ,s 3 ,si) 
1 



K 



(— 



H 



(— 



V2 



3(K(si, s 3 , s 2 ) - K(s 2 , s 3 , si)) 



+ H(s 1 ,s 3 ,s 2 ) - H(s 2 ,s 3 ,s 1 ) 

3. Numerical results for the structure functions 
are given in Figs. 1-3. We plot si, s 2 averaged 
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functions wx{Q 2 



= < 



/ dsids 2 W x {Q 2 ,si,s 2 ) 
for X = A, C 

J ds 1 ds 2 sign(si - s 2 ) W X (Q 2 , si, s 2 ) 
for X = D,E 



where 

Corrections to the tree level results soon be- 
come very important. Because of this, and be- 
cause of the obvious difficulty in accessing experi- 
mentally the region close to threshold, significant 
direct tests of these predictions in the low en- 
ergy region will require very high statistics. The 
0(p i ) chiral prediction for the branching ratio 

for (t -> v r + (3tt)(Q 2 < QLx)) is 01 the or- 
der of 10~ 5 for Q max = 600 MeV and 10~ 4 for 
Qmax = 700 MeV [g. This shows that experimen- 
tal tests in the low energy region are difficult but 
not hopeless. Moreover, some of the qualitative 
features of the structure functions, which the chi- 
ral expansion allows to understand in detail, will 
also have consequences at higher energies and so 
can be tested experimentally more easily. 

The first feature, which the chiral expansion 
predicts, is that w A + ^ 



( — +) ^ (oo-) 

W n W It I , 



,(00-) 



C ~ A 

c . In fact, the four plots can not be distin- 
guished from each other within the resolution of 
Fig. |l|. The equality of the structure functions 
for the two charge modes indicates that wa and 
wc are dominated by the [210] partition j|. We 
have explicitly verified that this is the case. The 
near equality of wa and wc at low energies can 
be understood algebraically. In Eq. ( |l2] ) we read 
that 



W A - W C = 8xj\H\ 2 



(15) 

H vanishes at tree level according to (|l^). Thus 
the difference starts as the square of a quantity 
of 0(p 2 ). Furthermore, H is kinematically sup- 
pressed at low energies, because it vanishes for 
si = S2- Therefore wa ~ wc should be valid 
even at energies well above those where one would 
trust a one loop CHPT calculation. 

In Fig. |l|we also compare to the prediction from 
the vector meson dominance model (VMD) in || . 



Figure 1. Integrated structure function wa(Q 2 )- 
CHPT prediction at tree level (dashed), one loop 
(solid) and from a vector meson dominance model 
(dotted). The four functions, wa and wc for both 
modes 27r~7r + and 2tt°tt^ all look identical within 
the resolution of this diagram. 
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Below \/Cp of about 600 MeV, where we would 
trust CHPT, we find very good agreement be- 
tween the VMD model and the one loop chiral 
prediction. 

4. In Fig. H and H we plot wd(Q 2 ) and we(Q 2 )- 
These are much smaller than wa , and correspond- 
ingly more difficult to be measured. A measure- 
ment of them in the low energy region can only 
be considered at a b or r-charm factory. 

The behaviour of these structure functions is 
rather interesting. The one loop predictions from 
CHPT differ strongly for the two charge modes. 
The origin of this sign difference can again be 
understood in detail: neglecting \H\ 2 that as dis- 
cussed above is tiny near threshold, Wd and We 
are the real and imaginary part of the same func- 
tion: 



4x 3 (x 1 +x 2 )KH* 



(16) 



The leading contribution to this function comes 
from the interference of the tree level part of K 
and the one loop part of H. At tree level K ex- 
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Figure 2. Integrated structure function wd(Q 2 )- 
CHPT at one loop for 2tt~tt + (solid) and for 
27r°7r~ (dashed-dotted) and from the VMD model 
in H (dotted, identical prediction for both charge 
modes) 
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Figure 3. Integrated structure function we(Q 2 )- 
CHPT at one loop for 2tt~tt + (solid) and for 
27r°7r~ (dashed-dotted) and from the VMD model 
in H (dotted, identical prediction for both charge 
modes) 
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hibits a simple sign difference for the two charge 
modes: 

*<-+> = ^(-2/3 + Ob 2 )) 

= !^(2/3 + 0( P 2 )). (17) 

Since H does not change sign between the two 
charge modes the same sign difference as in K 
at tree level is reproduced in the two structure 
functions Wd and We- In our numerical calcu- 
lation the sign difference remains up to \fCp- ~ 
0.6 GeV, where large corrections to this calcula- 
tion are already expected. 

Note that CHPT predicts w^~ +) < and 
w E + ^ > 0, which should be correct near thresh- 
old. However, experimental data are available for 
Q 2 > 0.8 GeV 2 f§, which indicate w ( D ~ +) > 

and w E + ' < 0. Thus, unless the higher or- 
ders in the chiral expansion completely change 
the CHPT result, somewhere between threshold 
and Q 2 ~ 0.8 GeV 2 there must be a zero for both 
structure functions in the 2ir~ir + mode. This fea- 
ture is absent in all the phenomenological models 
in the literature of which we are aware. We can 
not predict from the chiral expansion, whether 
these zeros will occur at very low energies or only 
a little below Q 2 = [M p + M„) 2 w 0.8 GeV 2 . We 
suggest a careful experimental examination of wd 
and we from Q 2 = 0.8 GeV 2 downwards. 

5. Let us summarize our main results. First, 
we predict that wa and wc are very similar at 
threshold. This has already been noted in phe- 
nomenological models, but we have shown that 
wa — wc is actually a low energy theorem of 
QCD, which receives corrections only at next-to- 
next-to leading order. So the symmetry structure 
of QCD requires wa ~ wq well beyond the very 
low energy region. 

Second, near threshold we find a sign differ- 
ence between the two charge modes for wd and 
we- Taking into account also the available exper- 
imental data in 2ir~ir + channel we conclude that 
in this channel both structure functions should 
have a zero somewhere between threshold and 
Q 2 ~ 0.8 GeV 2 . This feature is absent in any 
of the phenomenological models of which we are 
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aware, and its experimental verification would be 
the first evidence for the presence of the [300] par- 
tition state [U in tau decays into three pions. 

With the same method we can calculate also 
the scalar part of the matrix element close to 
threshold. It turns out that, compared to the one 
loop prediction from CHPT, this scalar part is 
underestimated, if present at all, in phenomeno- 
logical models. A detailed analysis of this and 
other topics, such as the application of CHPT to 
decays into two pions, is given elsewhere ||. 

4. Issues at Higher Q 2 : Resonances 

1. 7T~7r° final states. There are some dis- 
crepancies between experimental data and predic- 
tions using the "chirally normalized vector meson 
dominance" , or "CN-VMD" model for %Kw final 
states (t]||. Prompted by this, we have started 
to systematically test the assumptions of the CN- 
VMD model, beginning with the simplest mode, 
i.e. the two pion final state. 

A completely general ansatz for the relevant 
form factor Fv(Q 2 ) is 



F V (Q 2 ) 



f P (Q 2 )g P .AQ 2 ) 



Ml 



BW,(Q 2 ) 



fp'(Q 2 )g P >™(Q 2 ) 



M 2 -BW P ,(Q 2 ) 



(18) 



where the • • • may indicate both higher resonances 
(p" , . . . ) and possible non-resonant contribu- 
tions. We use Breit-Wigner resonance factors 
normalized to BWx(O) = 1, 



BW X (Q 2 



Ml 



Q 2 



(19) 



The on-shell values of the various couplings can 
(approximately) be determined from experimen- 
tal data, 



\9 P UM 2 )\ 
\9p>™(M})\ 
\f P (M 2 )\ 
\fAMp.)\ 



= 6.05 

= 1.39 

= 0.117 GcV 2 

= 0.18 GeV 2 



(20) 



If we assume that the meson couplings are ap- 
proximately constant from Q 2 — up to the rel- 
evant resonance masses, we find 



Fv(Q 2 ) 



f P (M 2 )g p ^(M 2 ) 



M 2 p 



BW P (Q 2 ) 



1.20 

fAM 2 p ,)g P ^{M 2 ) 
M 2 



BW P '(Q 2 ) 



-0.12 



(21) 



Here we have used the experimental knowledge 
about the negative relative phase between the p 
and the p' contributions. Note that no reliable 
error estimates are possible here. For the value 
f P '(M 2 )g p ,^(M 2 )/M 2 = 0.12 we get an error 
of ±0.02 from the uncertainty in M p i, but the 
uncertainty from gpi^n is unknown and might po- 
tentially be larger. 

We will now compare these results with the 
CN-VMD model. If we dominate with the p only, 
we get 

F v (Q 2 ) = lxBW p (Q 2 ) (22) 

Compared to (pil|), we find a 20% discrepancy in 
the strength of the p contribution. This has al- 
ready been observed in ||. If, however, we in- 
clude the first two resonances, we get 



F V (Q 2 



1 + /3 




xBW p {Q 2 ) 



xBW p ,{Q 2 ) 



(23) 



where we have taken (5 from the fit N = 1 in Q] , 
and with three resonances, we get (fit N = 2 in 
§) 



F V (Q 2 



1 



1 + /3 + 7 

1.16 

P 

1 + /3 + 7 



xBW p (Q 2 ) 



xBW p ,(Q 2 ) 



7 



+ 



7 



1 + /3 + 7 



-0.04 



xBW p »(Q 2 



(24) 



These numbers show that the CN-VMD model 
compares reasonably well with (fjl|), if the first 
two or three p like states are included. This gives 
us some confidence that the CN-VMD model is 
indeed a good approximation. A more detailed 
study is in preparation pi . 

2. Kir final states. There are two such final 
states, K ir~ and K~ir°, which are related by 
isospin. Due to m s > m u ,md, both Fy and Fs 
may contribute |Oj. Relevant resonances with 
the correct quantum numbers are: 

J P = \-{F V ): K*{m2) and IT* (1410) 
J p = 0+ (F s ) : 1^(1430) 

In jl2| , we have constructed a model with includes 
both the if*(1410) and the 1^(1430). 

The relative strength (3k> ~ —0.135 of the 
AT* (1410) in the vector channel was estimated 
from BK(Ktt) (note that (3k' ~ ftp, which is 
consistent with SU(3) flavor symmetry expec- 
tations), and the strength of the ^(1430) by 
matching to 0(p 4 ) Chiral Perturbation Theory. 
We found that the iT"(1410) and the i4T *(1430) 
contribute with comparable strength (both about 
5% to the rate). Thus one has to disentangle 
them using angular distributions, i.e. by measur- 
ing Wb 7 Wsa and Wsf- 

3. 37r final states. Some important tasks in 
these final states are to establish the presence of 
a scalar contribution and to measure its size, and 
to put upper limits on the p' contribution in two 
pion subresonances. 

4. KttK and ttKtt final states. As opposed to 
the case of the 37r, here both the vector and the 
axial vector currents can contribute (i.e. both the 
ai and the p', or both the K\ and the K*'). The- 
oretical predictions for their ratio V : A vary con- 
siderably |3||@@,||. Thus it should be mea- 
sured, using angular distributions. Understand- 
ing these decay modes may turn out to be very 
useful for measurements of possible non standard 
model CP-violation p6| . This is because suitable 
interference terms between the vector and the ax- 
ial vector contribution can be measured without 



polarized taus, and without knowledge of the r 
rest frame; and because the two interfering am- 
plitudes (vector and axial vector) may have com- 
parable strengths. 

Another interesting issue regarding the ttKtt 
final states are the two axial resonances with 
strangeness, the ifi(1270) and the i^i(1400), 
which both can contribute here. 

We believe that there are good reasons to sus- 
pect that the K\ widths given in the particle data 
book |l7j might be considerably too small ||. 
These reasons are: (i) Otherwise the BR(ttKtt) 
within the CN-VMD model turn out to be about 
a factor of 2 • • • 3 too large. Of course, the obvi- 
ous thought is that this might be due to a failure 
of the CN-VMD model. However, due to the fact 
that most of the other predictions of the model 
for three meson final states agree reasonably well 
with data, and due to the various symmetry re- 
lations between the different three meson final 
states, we are actually unable to find any other 
natural explanation. This led us to suspect that 
indeed the true K\ widths might be larger, (ii) 
Given that T ai « 400 • • • 600 MeV, the K x widths 
in p7| appear unnaturally small, (iii) Experi- 
ments on which the numbers quoted in JlTj are 
based, have made the same assumptions for the 
parametrization of background amplitudes, which 
in the past yielded very small values for T ai , and 
which subsequently were superseded by cleaner 
measurements of r ai using r decays. Thus our 
conclusion is that one should measure the K\ res- 
onance parameters in r — * ttKhv t . In Fig. ^ 
we show that indeed u>b is very sensitive to the 
K\ parameters. These issues will be discussed in 
more detail in a forthcoming publication poj . 

5. Conclusions 

Let us conclude by making some suggestions 
to experimentalists about what we believe to be 
interesting and feasible measurements. 

First, let us consider decays with three pions. 
We would suggest to measure wd and we for both 
charge modes, 2tt^tt + and 2tt°tt~. These two 
structure functions might very well be different 
for the two modes. This would be the first evi- 
dence for the [300] partition class, which decays 
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Figure 4. Structure functions for r — » X 7r 7r + : 
wa(Q 2 ) (solid) and w s (Q 2 ) (dashed). The 
larger solid curve with two clearly distinguishable 
peaks is obtained using the RPP-96 values for 
IVi(i27o) an d IV^uoo), the smaller solid curve 
using rV l(12 7o) = r Kl ( 1400 ) = 250 MeV 
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differently into 2tt~it + and 27r°7r~. Furthermore, 
we would like to urge to measure Wd and we for 
Q 2 less or nearly equal to 0.8 GeV 2 for the 2tt~7t + 
charge mode. There is most probably a zero and 
a change of sign both in wd and we in this mode, 
somewhere between Q 2 = 9M 2 ■ • • 0.8 GeV 2 . 

We also suggest to measure the scalar contri- 
bution via wsb, and to try to put limits on the 
p' contribution in the two-pion subresonances. 

Secondly, regarding r — > irKirv T , we suggest 
to measure the parameters of the K\ resonances 
using this r decay mode. 
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